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Abstract A search for supersymmetry (SUSY) in events 
with large missing transverse momentum, jets, and at 
least one hadronically decaying r lepton, with zero or 
one additional light lepton (e//i), has been performed 
using 4.7 fb~^ of proton-proton collision data at ^/s — 
7 TeV recorded with the ATLAS detector at the Large 
Hadron Collider. No excess above the Standard Model 
background expectation is observed and a 95% con- 
fidence level visible cross-section upper limit for new 
phenomena is set. In the framework of gauge-mediated 
SUSY-breaking models, lower limits on the mass scale 
A are set at 54 TeV in the regions where the fi is the 
next-to- lightest SUSY particle (tan (3 > 20). These lim- 
its provide the most stringent tests to date of GMSB 
models in a large part of the parameter space consid- 
ered. 
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1 Introduction 

This paper reports on the search for supersymmetry 
(SUSY) [IHS] in events with large missing transverse 
momentum, jets and at least one hadronically decaying 
r lepton. Four different topologies with a r in the final 
state have been studied: one r lepton, at least two r 
leptons, one t lepton and precisely one additional muon 
and one r lepton and precisely one additional electron. 
The minimal gauge-mediated supersymmetry-breaking 
(GMSB) model [TUUTF is considered as benchmark to 
evaluate the reach of this analysis. 

SUSY introduces a symmetry between fermions and 
bosons, resulting in a SUSY partner (sparticle) for each 
Standard Model (SM) particle with identical mass and 
quantum numbers except a difference by half a unit 



of spin. Assuming i?-parity conservation [T6Vl20j . spar- 
ticles are produced in pairs. These would then decay 
through cascades involving other sparticles until the 
lightest SUSY particle (LSP), which is stable, is pro- 
duced. Since equal mass SUSY partners are excluded, 
SUSY must be a broken symmetry. Minimal GMSB 
models can be described by six parameters: the SUSY- 
breaking mass scale in the low-energy sector {A), the 
messenger mass (Mmoss), the number of SU(5) messen- 
ger fields {N5), the ratio of the vacuum expectation val- 
ues of the two Higgs doublets (tan /3) , the Higgs-sector 
mixing parameter (fi) and the scale factor for the grav- 
itino mass (Cgrav)- Foi' the analysis presented in this pa- 
per, A and tan /? are treated as free parameters, and the 
other parameters are fixed to the values already used 
in Refs. [211 Hi]: M^oss = 250 TeV, A5 = 3, ^ > and 
Cgrav = 1- The Cgrav parameter determines the lifetime 
of next-to-hghtest SUSY particle (NLSP); for Cgrav = 1 
the NLSP decays promptly (ctnlsp < 0.1mm). With 
this choice of parameters, at moderate A the production 
of gluino and/or squark pairs is expected to dominate 
at the LHC; these sparticles will decay into the next-to- 
lightest SUSY particle (NLSP), which subsequently de- 
cays to the LSP. In GMSB models, the LSP is the very 
light gravitino (G). The NLSP is the dominant sparti- 
cle decaying to the LSP and this leads to experimental 
signatures which are largely determined by the nature 
of the NLSP. This can be either the lightest stau (fi ) , a 
right-handed slepton {£r), the lightest neutralino (xi), 
or a sneutrino (;/), dominantly leading to final states 
containing r leptons, light leptons {£ = e,/x), photons, 
6-jets, or neutrinos. At large values of tan/3, the fi is 
the NLSP for most of the parameter space, which leads 
to final states containing at least two r leptons. In the 
so-called CoNLSP region, where the mass difference be- 
tween the fi and the In is smaller than the sum of the 
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T and light-lepton masses, both the fi and the £_r decay 
directly into the LSP and are therefore NLSPs. 

Previous searches for fi pair production, with the 
subsequent decay f i — > tG in the minimal GMSB model, 
have been reported by the LEP Collaborations ALEPH 
[23], DELPHI [24 and OPAL [25,. The analysis re- 
ported in this paper extends the searches in 2fb^^ of 
data presented in Refs. [2T1 [32]. It comprises the full 
2011 dataset, corresponding to an integrated luminosity 
of (4.7 ± 0.1) fb-i [261 [22 after applying beam, detec- 
tor and data-quality requirements. A complementary 
search interpreted in GMSB, requiring two light lep- 
tons, has also been performed using the same dataset 
by the ATLAS coUaboration [25] ■ The CMS Collabora- 
tion has searched for new phenomena in same-sign r- 
pair events | 29j and multi-lepton events including two 
r leptons in the final state [3D] using 35pb~^ of data, 
where the minimal GMSB model was not considered. 



2 ATLAS detector 

The ATLAS experiment [31] is a multi-purpose detec- 
tor with a forward-backward symmetric cylindrical ge- 
ometry and nearly 47r solid angle coverage. The inner 
tracking detector (ID) consists of a silicon pixel detec- 
tor, a silicon microstrip detector and a transition radi- 
ation tracker. The ID is surrounded by a thin super- 
conducting solenoid providing a 2 T magnetic field and 
by fine-granularity lead/liquid-argon (LAr) electromag- 
netic calorimeters. An iron/scintillator-tile calorimeter 
provides hadronic coverage in the central pseudorapid- 
it}Q range. The endcap and forward regions are instru- 
mented with liquid-argon calorimeters for both elec- 
tromagnetic and hadronic measurements. An extensive 
muon spectrometer system that incorporates large su- 
perconducting toroidal magnets surrounds the calorime- 
ters. 



3 Simulated samples 

The Monte Carlo (MC) simulations used to evaluate 
the expected backgrounds and selection efficiencies for 
the SUSY models considered are very similar to the 
ones used in Refs. [H] [22]. A suite of generators is 

^ ATLAS uses a right-handed coordinate system with its ori- 
gin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The a;-axis 
points from the IP to the centre of the LHC ring and the 
2/-axis points upward. Cylindrical coordinates (r, (f>) are used 
in the transverse plane, (f> being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 6 as r/ = — lntan(0/2). 



used to aid in the estimate of SM background contri- 
butions. The ALPGEN generator [32] is used to simulate 
samples of W and ^/7* events with up to five (for Z 
events) or six (for W events) accompanying jets, where 
CTEQ6L1 [33^ is used for the parton distribution func- 
tions (PDFs). Z/j* events with mu < 40 GeV are re- 
ferred to in this paper as "Drell-Yan" . Top quark pair 
production, single top production and diboson {WW 
and WZ) pair production are simulated with MCONLO 
[34H36] and the next-to-leading-order (NLO) PDF set 
CTIO [37]. Fragmentation and hadronization are per- 
formed with Herwig ^5], using JIMMY [32] for the un- 
derlying event simulation. The decay of r leptons and 
radiation of photons are simulated using TAUOLA [TOIBT] 
and PHOTOS [42j , respectively. The production of multi- 
jet events is simulated with PYTHIA 6.4.25 [13] using 
the AUET2B tune 44J and MRST2007 LO* ^ PDFs. 
The SUSY mass spectra are calculated using ISAJET 
7.80 [46]. The MC signal samples are produced using 
Herwig++ 2.4.2 47J with MRST2007 LD* PDFs. Signal 
cross-sections are calculated to next-to-leading order in 
the strong coupling constant, adding the resummation 
of soft gluon emission at next-to-leading-logarithmic ac- 
curacy (NLO-hNLL) [15H52]- The nominal SUSY pro- 
duction cross-sections and their uncertainties are taken 
from an envelope of cross-section predictions using dif- 
ferent PDF sets and factorisation and renormalization 
scales, as described in Ref. [53]. The GMSB signal sam- 
ples are generated on a grid ranging from yl = 10 TeV 
to yl = 80 TeV and from tan /3 = 2 to tan /3 = 67, with 
the cross-section dropping from 100 pb for A — 15 TeV 
to 5.0 fb for vl = 80 TeV. 

All samples are processed through the Geant4- 
based simulation [SJ of the ATLAS detector [55j . The 
full simulation also includes a realistic treatment of the 
variation of the number of pp interactions per bunch 
crossing (pile-up) in the data, with an average of nine 
interactions per crossing. 

4 Object reconstruction 

Jets are reconstructed using the anti-fct jet clustering 
algorithm [56j with radius parameter R = 0.4. Jet en- 
ergies are calibrated to correct for upstream material, 
calorimeter non-compensation, pile-up, and other ef- 
fects [57 . Jets are required to have transverse momenta 
(pt) greater than 25 GeV and |?7| < 2.8, except in the 
computation of the missing transverse momentum, where 
|?7| < 4.5 and pt greater than 20 GeV is required. 

Muon candidates are identified as tracks in the ID 
matched to track segments in the muon spectrometer 
[55] . They are required to have pr > 10 GeV and \tj\ < 
2.4. Electron candidates are constructed by matching 
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Table 1 Event selection for the four final states presented in this paper. Numbers in parentheses are the minimum transverse 
momenta required for the objects. Pairs of numbers separated by a slash denote different selection criteria imposed in different 
data-taking periods. 



It I 2t I T-I-/1 I T+e 



Trigger 


jetMET 
pi^* > 75 GeV 
E;^''''' > 45/55 GeV 


jetMET 
pi^' > 75 GeV 
E;^'""" > 45/55 GeV 


muon /muon+jet 
p!^ > 18 GeV 
pj^' > 10 GeV 


electron 

> 20/22 GeV 


Jet req. 
i;^?'^^ req. 


>2jets (130, 30 GeV) 
E^^'"^" > 130/150 GeV 


= 1 medium (20 GeV), 
=0 loose 


>2jets (130, 30 GeV) 
> 130/150 GeV 


> 2 loose (20 GeV) 


>ljet (50 GeV) 

1 (20 GeV) 
>1 medium (20 GeV) 


1 e (25 GeV) 
>1 medium (20 GeV) 


Kinematic 
criteria 


^(0j.t,,-p»"-.) >0.3; 
E^^^^ym^fi > 0.3, 
rriT > 110 GeV; 
Ht > 775 GeV 


^(</'jet,,-p™-) >0.3 

m'^ + mJ^ > 100 GeV 
Ht > 650 GeV 


> 100 GeV 
meff > 1000 GeV 


^ > 100 GeV 
meff > 1000 GeV 



electromagnetic clusters with tracks in the ID. They 
are then required to satisfy px > 20 GeV, \r]\ < 2.47 
and to pass the "tight" identification criteria described 
in Ref. [59], re-optimized for 2011 conditions. 

Electrons or muons are required to be isolated, i.e. 
the scalar sum of the transverse momenta of tracks 
within a cone of AR = ^(zi</«)2 -I- {Arjy < 0.2 around 
the lepton candidate, excluding the lepton candidate 
track itself, must be less than 10% of the lepton's trans- 
verse energy for electrons and less than 1.8 GeV for 
muons. Tracks selected for the electron and muon isola- 
tion requirement defined above have px > 1 GeV and 
are associated to the primary vertex of the event. 

The missing transverse momentum vector p^^^^ (and 
its magnitude iJ™"'*') is measured from the transverse 
momenta of identified jets, electrons, muons and all 
calorimeter clusters with |?7| < 4.5 not associated to 
such objects [BO.- For the purpose of the measurement 
of Et^^^, t leptons are not distinguished from jets. 

Jets originating from decays of 6-quarks are iden- 
tified and used to separate the W and tt background 
contributions. They are identified by a neural-network- 
based algorithm, which combines information from the 
track impact parameters with a search for decay ver- 
tices along the jet axis [GT] . A working point corre- 
sponding to 60 % tagging efficiency for 6-jets and < 1 % 
mis-identification of light-flavour or gluon jets is chosen 
[62]. 

The r leptons considered in this search are recon- 
structed through their hadronic decays. The r recon- 
struction is seeded from anti-fct jets {R = 0.4) with 
Pt > 10 GeV. An rj- and px-dependent energy calibra- 
tion to the hadronic r energy scale is applied. Discrimi- 
nating variables based on track information and observ- 



ables sensitive to the transverse and longitudinal shape 
of the energy deposits of r candidates in the calorimeter 
are used. These quantities are combined in a boosted 
decision tree (BDT) discriminator [63] to optimize their 
impact. Calorimeter information and measurements of 
transition radiation are used to veto electrons mis-iden- 
tified as T leptons. Suitable r lepton candidates must 
satisfy pr > 20 GeV, |?7| < 2.5, and have one or three 
associated tracks of px > 1 GeV with a charge sum of 
±1. A sample of Z —> rr events is used to measure the 
efficiency of the BDT r identification. The "loose" and 
"medium" working points in Ref. |63j are used herein 
and correspond to efficiencies of about 60 % and 40 % 
respectively, independent of pxj with a rejection factor 
of 20-50 against r candidates built from hadronic jets 
("fake" r leptons). 

5 Event Selection 

Four mutually exclusive final states are considered for 
this search: events with only one "medium" r, no addi- 
tional "loose" T candidates and no muons or electrons, 
referred to as 'Ir'; events with two or more "loose" r 
candidates and no muons or electrons, referred to as 
'2r'; events with at least one "medium" r and exactly 
one muon ('t+^') or electron ('r-fe'). 

In the It and 2t final states, candidate events are 
triggered by requiring a jet with high transverse mo- 
mentum and high S^''^" ('jetMET') [65], both measured 
at the electromagnetic scal^ In the r+/i final state, 

^The electromagnetic scale is the basic calorimeter signal 
scale for the ATLAS calorimeters. It has been established 
using test-beam measurements for electrons and muons to 
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Table 2 Definition of the background control regions (CRs) used to estimate the normalization of background samples in the 
four final states: It, 2t, r+// and r+e. 



Background | It | 2t | t+/i | r + e 



tt 


^(<7ijoti.,-p-'-) > 0.3rad 
niT < 70 GeV 
E^''"'/m^ii > 0.3 
6-tag template fit 


^(■Jf-jct^.^-p"-) >0.3rad 
m'^^ + ml^ > 100 GeV 
Ht < 550 GeV 

A^b-tag > 1 


SO GeV < E^'"^" < 100 GeV 
50 GeV < ^ < 150 GeV 

A^b-tag > 1 


W+iets 


^(<j!'jeti,-p»i-) > O.Srad 
TUT < 70 GeV 
B^f'^Vmeff > 0.3 


^(■^jet,,-p»i™) > 0.3rad 
m'Jj} + m:p > 100 GeV 
Ht < 550 GeV 

Affc-tag = 


SO GeV < < 100 GeV 
50 GeV < m^'^ < 150 GeV 
Ar^_tag = 


Z+jets 


2/i (20 GeV), |r?| < 2.4 
>2jets (130, 30 GeV) 

Afi,_tag = 


^('?ijet,,,-p--0 > O.Srad 
m:^' + m^^ < 80 GeV 
i?T < 550 GeV 


MC-based normalization 


Multi-jet 


^('^jeti,-p»i-) < 0.3rad 
_EJf,"'=Vmeff < 0.3 


^{<t>iot,^-pg^'') < O.Srad 
£™'=7moff < 0.4 


compare events with and without 
lepton isolation [64] 



events are selected by a muon trigger and a muon-plus- 
jet trigger ('muon+jet'), while in the r+e final state, a 
single-electron trigger requirement is imposed [6 5) . The 
trigger requirements have been optimized to ensure a 
uniform trigger efiiciency for all data-taking periods, 
which exceeds 98 % with respect to the offline selection 
for all final states considered. 

Pre-selected events are required to have a recon- 
structed primary vertex with at least five tracks (with 
Pt > 0.4 GeV). To suppress soft multi-jet events in the 
It and 2r final states, a second jet with > 30 GeV is 
required. Remaining multi-jet events, where highly en- 
ergetic jets are mis-measured, are suppressed by requir- 
ing the azimuthal angle between the missing transverse 
momentum vector and either of the two leading jets to 
be greater than 0.3 rad. Three quantities characterising 
the kinematic properties of the event are used to fur- 
ther suppress the main background processes (VF-|-jets, 
Z-l-jets and tt events) in all four final states: 

— the transverse mass formed by E^^'^^ and either 
the pt of the t lepton in the Ir and 2r channels, or 
of the light lepton (e//i) in the r+fj, and T+e ones: 

m^'^ = ^J2p''/E^''^%l - cos{A(j){T/£, £;™""))) ; 

— the scalar sum _ffx of the transverse momenta of r 
lepton candidates and the two highest momentum 
jets in the events: Ht = J2Pt + Si^i ' ' 

— the effective mass rrics = Ht + Et^^^- 



give the correct response for the energy deposited in electro- 
magnetic showers, although it does not correct for the lower 
response of the calorimeter to hadrons. 



For each of the four final states, specific criteria are ap- 
plied to the above quantities in order to define a signal 
region (SR), as summarized in Table [l] 

Figure [l] shows the and mJ^ distributions 
for the It and 2r channels after all the requirements of 
the analysis except the final requirement on Ht- Sim- 
ilarly, Figure [2] shows the distributions for the 
T+fj, and T+e channels after all the requirements of 
the analysis except the final mes requirement. 

Figure [3] and |4] show the Ht distributions in the 
It and 2t channels, and mog distributions in the t+h 
and T+e channels, respectively, after all other selection 
criteria have been imposed. 

6 Background estimation 

The SM background expectation predicted by simula- 
tion in the SR is corrected by means of control regions 
(CRs), which are chosen such that a specific background 
process is enriched while any overlap with the SR is 
avoided. Data/MC comparison in the CRs show that 
MC overestimates the number of events compared to 
data, mainly due to mis-modelling of r mis-identifica- 
tion probabilities and kinematics. Scaling factors are 
therefore obtained from the ratio of the number of ob- 
served events to the number of simulated background 
events in the control region where a given background 
contribution is enriched. Studies comparing data with 
MC simulations show that the t mis-identification prob- 
ability is, to a good approximation, independent of the 
kinematic variables used to separate the SR from the 
CRs, so that the measured ratio of the data to MC 
event yields in the CR can be used to compute scaling 
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(a) rriT distribution for the Ir final state. 
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Fig. 1 Distribution of (a) mx and (b) raJ^ + ml^ for the 
It and 2t final states, respectively, after all analysis require- 
ments but the final requirement on Ht ■ Data are represented 
by the points, with statistical uncertainty only. The SM pre- 
diction includes the data-driven corrections discussed in the 
text. The band centred around the total SM background in- 
dicates the uncertainty due to finite MC sample sizes on the 
background expectation. Also shown is the expected signal 
from two typical GMSB samples {A = 50 TeV, tan /3 = 40, 
A = 50 TeV, tan /3 = 20). 



factors to correct the MC background prediction in the 
SR. 

The dominant background contributions in the SR 
arise from top quark pair and single top events (here- 
after generically indicated as 'top'), VF+jets, Z+jets 
and multi-jet events. The latter background does not 
contribute significantly to the t+h final state. The CR 
definitions used to estimate these background contribu- 
tions in the various channels are summarized in Table 

H 



Fig. 2 Distribution of mf^ ^ for the (a) r-|-/i and (b) T + e fi- 
nal states after all analysis requirements but the final require- 
ment on rrieff . Data are represented by the points, with statis- 
tical uncertainty only. The SM prediction includes the data- 
driven corrections discussed in the text. The band centred 
around the total SM background indicates the uncertainty 
due to finite MC sample sizes on the background expectation. 
Also shown is the expected signal from two typical GMSB 
samples {A = 50 TeV, tan /3 = 40, A = 50 TeV, tan /3 = 20). 



6.1 Background estimation in the 2r channel 

The W and top background contributions are domi- 
nated by events in which one r candidate is a true r and 
the others are mis-reconstructed from hadronic activity 
in the final state. The background from Z+jets events is 
dominated by final states with Z — > tt decays. The CRs 
defined for the estimation of these background contri- 
butions have a very small contamination from multi-jet 
events due to the requirement on Z\((/)jj,(.^ ^_pmiss) and 
the presence of two or more r leptons. The signal contri- 
bution in these CRs is expected to be at less than 0.1 % 
for the models considered. Correlations between differ- 
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Fig. 3 Distribution of Ht for the (a) It and (b) 2t final 
states after all analysis requirements. Data are represented 
by the points, with statistical uncertainty only. The SM pre- 
diction includes the data-driven corrections discussed in the 
text. The band centred around the total SM background in- 
dicates the uncertainty due to finite MC sample sizes on the 
background expectation. Also shown is the expected signal 
from two typical GMSB samples {A = 50 TeV, tan /3 = 40, 
A = 50 TeV, tan /3 = 20). 

ent samples in the various CRs are taken into account 
by considering the matrix equation N'^''*'* — A us, 
where N'^^'^ is the observed number of data events in 
each of the CRs defined in Table[2] after subtracting the 
expected number of multi-jet events and any remaining 
sub-dominant background contribution, obtained from 
MC simulation. The matrix A is obtained from the MC 
expectation for the number of events originating from 
each of the background contributions (top, W and Z). 
The vector a; of scaling factors is then computed by 
inverting the matrix A. To obtain the uncertainties for 
the scaling factors, all contributing parameters are var- 
ied according to their uncertainties, the procedure is re- 



Fig. 4 Distribution of meff for the (a) r-h/x and (b) r + e final 
states after all analysis requirements. Data are represented 
by the points, with statistical uncertainty only. The SM pre- 
diction includes the data-driven corrections discussed in the 
text. The band centred around the total SM background in- 
dicates the uncertainty due to finite MC sample sizes on the 
background expectation. Also shown is the expected signal 
from two typical GMSB samples {A = 50 TeV, tan /3 = 40, 
yl = 50 TeV, tan /3 = 20). In the top figure, the event in data 
surviving all the analysis requirements is shown in the over- 
flow bin. 



peated and new scaling factors are obtained. The width 
of the distribution of each resulting scaling factor is 
used as its uncertainty. The typical scaling factors ob- 
tained with this procedure are between 0.75 and 1 , with 
uncertainty of order 40 %. The multi-jet background ex- 
pectation is computed in a multi-jet-dominated CR de- 
fined by inverting the Zi((/)j(,j^ ^_pmias) requirement and 
not applying the m^^ -I- ml^ and Ht selection. In addi- 
tion, an upper limit is imposed on the ratio E^^^^^/nieS 
to increase the purity of this CR sample. 
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6.2 Background estimation in the It channel 

The number of events from VF+jets and WZ processes 
in the SR is estimated by scahng the number of cor- 
responding MC events with the ratio of data to MC 
events in the VF+jets CR. The corresponding scahng 
factors are computed separately for the cases in which 
the r candidates from Vt^/top decays are true r lep- 
tons and for those in which they are mis-reconstructed 
from hadronic activity in the final state. It has been 
checked that the same scaling factors can be applied to 
both M^-l-jets and WZ processes. In the case of W^-l-jets 
background events with true r candidates, the charge 
asymmetry method [B5J [HZ] is used. To estimate the 
background from top events with true r candidates, a 
scaling-factor-based-technique is also used, where the 
number of 6-tagged events in data in the top CR is fit- 
ted to a template from MC simulation ('template fit'). 
For background events in both W/top processes due 
to fake r candidates, the matrix method already dis- 
cussed for the 2r background estimation is employed, 
where the parameters in the vector a; of scaling factors 
are w^y^®, uJ^^ , '*^top'' ^^^d wJop*^. The region dominated 
by fake r candidates is defined by mx > 110 GeV and 
iJx < 600 GeV, while the one dominated by true r 
candidates is defined by requiring mx < 70 GeV. The 
values of w^™ obtained from this method and from 
the template fit are in very good agreement. The fac- 
tor obtained with the charge asymmetry method 
agrees within 2ct with the one obtained with the ma- 
trix inversion method. The difference between the two 
Wj^"'' values is then assigned as a systematic uncertainty 
on the W^-fjets background estimation procedure. The 
background from Z+jets events is due to events where 
the Z decays to a pair of neutrinos, and contributes 
fully to the observed E^^^^. The background contribu- 
tion in the SR is estimated from data by measuring the 
data/MC ratio from Z — )■ £~^£~ decays in the Z-|-jets 
CR defined in Table [2] Typical scaling factors are be- 
tween 0.75 and 1.2 , with uncertainty of order 20 %. The 
multi-jet background expectation is computed in the 
same way as in the 2t channel. 

6.3 Background estimation in the r+fi and r+e 
channels 

The top background contribution consists of events where 
the muon (electron) candidate is a true muon (elec- 
tron), and the r candidate can either be a true r or a 
hadronic jet mis-identified as a r. On the other hand, 
the W-|-jets background consists mainly of events where 
the T candidate is mis-reconstructed from hadronic ac- 
tivity in the final state. For this reason, the top CR 



is divided into two subregions: one dominated by true 
r candidates, defined by 100 GeV < m^''^ < 150 GeV, 
and one dominated by fake ones (50 GeV < ^ < 
100 GcV). The same matrix approach already described 
is then used to estimate the true/fake top and VF-|-jets 
background contributions to the SR. The scaling factors 
obtained are about 0.6-0.8 , with typical uncertainty of 
15 %. The Z-|-jets background is much smaller than the 
W^-l-jets one, and it is estimated using MC simulated 
events. The multi-jet background arises from mis-iden- 
tified prompt leptons. By comparing the rates of events 
with and without the lepton isolation requirement, a 
data-driven estimate is obtained following the method 
described in Ref . |M] . 

The contribution from other sources of background 
considered (Drell-Yan and diboson events) is estimated 
in all analyses using directly the MC normalizations, 
without applying any further scaling factor. 

Table [S] summarizes the estimated numbers of back- 
ground events in the SR for each channel. 



7 Systematic uncertainties on the background 

Various systematic uncertainties were studied and the 
effect on the number of expected background events 
in each channel presented was evaluated, following the 
approach of Refs. [2T1 [32]. The dominant systematic 
uncertainties in the different channels are summarized 
in Table m 

The theoretical uncertainty on the MC-based cor- 
rected extrapolation of the VF+jets and top backgrounds 
from the CR into the SR is estimated using alternative 
MC samples. These MC samples were obtained by vary- 
ing the renormalization and factorisation scales, the 
functional form of the factorisation scale and the match- 
ing threshold in the parton shower process in the gen- 
erators used for the simulation of the events described 
in section [3l 

Systematic uncertainties on the jet energy scale ( JES) 
and jet energy resolution ( JER) [57j are applied in MC 
events to the selected jets and propagated throughout 
the analysis. The difference in the number of expected 
background events obtained with the nominal MC sim- 
ulation after applying these changes is taken as the sys- 
tematic uncertainty. 

The effect of the r energy scale (TES) uncertainty 
on the expected background is estimated in a similar 
way. The uncertainties from the jet and r energy scale 
are treated as fully correlated. 

The uncertainties on the background estimation due 
to the r identification efficiency depend on the r identi- 
fication algorithm ( "loose" or "medium" ) , the kinemat- 



8 



Table 3 Number of expected background events and data yields in the four final states discussed. Where possible, the 
uncertainties are separated into statistical and systematic parts. The SM prediction is computed taking into account correlations 
between the different uncertainties. Also shown are the number of expected signal MC events for one GMSB point (yl=50 TeV, 
tan/3=20), the 95% confidence level (CL) upper limit on the number of observed (expected) signal events and corresponding 
cross-section from any new physics scenario that can be set for each of the four final states, taking into account the observed 
events in the data and the background expectations. 





Ir 


2r 


r + fi 


T + e 


Multi-jet 
W + jets 
Z + jets 
Top 

Diboson 
Drell-Yan 


0.17 ±0.04 ±0.11 
0.31 ±0.16 ±0.16 
0.22 ±0.22 ±0.09 
0.61 ±0.25 ±0.11 

< 0.05 

< 0.36 


0.17 ±0.15 ±0.36 
1.11 ±0.67 ±0.30 
0.36 ±0.26 ±0.35 
0.76 ±0.31 ±0.31 
0.02 ±0.01 ±0.07 
0.49 ±0.49 ±0.21 


< 0.01 
0.27 ±0.21 ±0.13 
0.05 ±0.05 ±0.01 
0.36 ±0.18 ±0.26 
0.11 ±0.04 ±0.02 

< 0.002 


0.22 ±0.30 
0.24 ±0.17 ±0.27 
0.17 ±0.12 ±0.05 
1.41 ±0.27 ±0.84 
0.26 ±0.12 ±0.11 
< 0.002 


Total background 


1.31 ±0.37 ±0.65 


2.91 ±0.89 ±0.76 


0.79 ±0.28 ±0.39 


2.31 ±0.40 ±1.40 


Signal MC Events 

{A = 50 TeV, tan /3 = 20) 


2.36 ±0.30 ±0.60 


4.94 ±0.45 ±0.74 


2.48 ±0.30 ±0.39 


4.21 ±0.38 ±0.46 


Data 


4 


1 


1 


3 




Obs. (exp.) upper limit on 
number of signal events 


7.7 (4.5) 


3.2 (4.7) 


3.7 (3.4) 


5.2 (4.6) 


Obs. (exp.) upper limit on 
visible Cross-Section (fb) 


1.67 (0.95) 


0.68 (0.99) 


0.78 (0.72) 


1.10 (0.98) 



ics of the T sample and the number of associated tracks. 
In the different channels, they vary between 2-5%. 

A systematic uncertainty associated with the sim- 
ulation of pile-up in the MC events is also taken into 
account, with uncertainties varying between 5-20%. 

The effect of the 1.8% uncertainty on the luminos- 
ity measurement [5^ [27] is also considered on the nor- 
malization of the background contributions for which 
scale factors derived from CR regions were not applied 
(Drell-Yan and diboson in all channels, and Z+jets in 
the T+/i and r + e channels). 

The total systematic uncertainties obtained in the 
It, 2r, r+/i and r+e channels are 52 %, 26 %, 49 % and 
60%, respectively. The limited size of the MC samples 
used for background estimation gives rise to a statistical 
error ranging from 21% in the It channel to 46% in 
the T+e channel. 



states, in the fi NLSP region the efficiency is about 3 % 
for the 2t channel, 1 % for the r+fj, and T + e channels, 
and 0.5% for the It channel. In the non-fi NLSP re- 
gions and for high A values it drops to 0.1-0.2% for 
all final states. The total systematic uncertainty on the 
signal selection from the various sources discussed in 
Section [7] ranges between 10-15 % for the It channel, 
15-18% for the 2t channel, 8-16% for the T+/i chan- 
nel and 11-17% for the T+e channel over the GMSB 
signal grid. 

Theoretical uncertainties related to the GMSB cross- 
section predictions are obtained using the same proce- 
dure as detailed in Ref. [22] • These uncertainties are cal- 
culated for individual SUSY production processes and 
for each model point in the GMSB grid, leading to over- 
all theoretical cross-section uncertainties between 5 % 
and 25%. 



8 Signal efficiencies and systematic 
uncertainties 

The GMSB signal samples are described in Section [3| 
The total cross-section drops from 100 pb for yl = 15 TeV 
to 5.0 fb for A = 80 TeV. The cross-section for strong 
production, for which this analysis has the largest effi- 
ciency, decreases faster than the cross-sections for slep- 
ton and gaugino production, such that for large val- 
ues of A the selection efficiency with respect to the to- 
tal SUSY production decreases. For the different final 



9 Results 

Table [3] summarizes the number of observed data events 
and the number of expected background events in the 
four channels, with separate statistical and systematic 
uncertainties. No significant excess is observed in any 
of the four signal regions. From the numbers of ob- 
served data events and expected background events, up- 
per limits at 95% confidence level (CL) of 7.7, 3.2, 3.7 
and 5.2 signal events from any scenario of physics be- 
yond the SM are calculated in the It, 2t, t+/x and T + e 
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Table 4 Overview of the major systematic uncertainties and the MC statistical uncertainty for the background estimates in 
the four channels presented in this paper. 



Source of Uncertainty 


Ir 


2t 




r+e 


CR to SR Extrapolation 


27% 


12% 


26% 


29% 


Jet Energy Resolution 


21% 


6.5% 


5.4% 


13% 


Jet Energy Scale 


20% 


4.8% 


11% 


8.5% 


T Energy Scale 


10% 


8.5% 


0.3% 


4.3% 


Pile-up modelling 


5.1% 


14% 


20% 


3.5% 


MC statistics 


21% 


32% 


39% 


46% 



channels, respectively. Using only the background pre- 
dictions, expected limits of 4.5, 4.7, 3.4 and 4.6 events 
are obtained for the four channels (Ir, 2r, r+fi and 
T+e). The limits on the number of signal events are 
computed using the profile likelihood method [SS] and 
the CLs criterion [51] . Uncertainties on the background 
and signal expectations are treated as Gaussian-dis- 
tributed nuisance parameters in the likelihood fit. The 
signal-event upper limits translate into a 95 % CL ob- 
served (expected) upper limit on the visible cross-sec- 
tion for new phenomena for each of the four final states, 
defined by the product of cross-section, branching frac- 
tion, acceptance and efficiency for the selections defined 
in Section |5] The results are summarized in Table [3] for 
all channels. In order to produce the strongest possi- 
ble 95 % CL limit on the GMSB model parameters A 
and tan/?, a statistical combination of the four chan- 
nels is performed. The likelihood function representing 
the outcome of the combination includes the statisti- 
cal independence of the four final states considered. 
The resulting observed and expected lower limits for 
the combination of the four final states are shown in 
Figure [5j These limits are calculated including all ex- 
perimental and theoretical uncertainties on the back- 
ground and signal expectations. Excluding the theoret- 
ical uncertainties on the signal cross-section from the 
limit calculation has a negligible effect on the limits 
obtained. Figure [5] also includes the limits from OPAL 
[25] for comparison. The best exclusion from the com- 
bination of all final states is obtained for yl = 58 TeV 
for values of tan (3 between 45 and 55 . The results ex- 
tend previous limits and values of yl < 54 TeV are now 
excluded at 95 % CL, in the regions where the fi is the 
next-to-hghtest SUSY particle (tan/3 > 20). 

10 Conclusions 

A search for SUSY in final states with jets, EJ^'^^, fight 
leptons (e//z) and hadronically decaying r leptons is 
performed using 4.7 fb~^ of ^/s = 7 TeV pp colhsion 
data recorded with the ATLAS detector at the LHC. 



GMSB: M„,ess=250 TeV, N5=3, ^>0, Cg„„=1 




f ? 
a a 

S 1 



^^^^H Observ 


ed 95% CL limit (±1 o^^^) 


Expected 95% CL limit (±1 a„p) 


2.05tb 


2t: ATLAS 95%CLiimit 


OPAL 


5%CL(i,) 




OPAL 


5% CL ( fi^j 




OPAL 


5% CL 





CoNLSP 



10 20 30 40 50 60 70 80 90 100 

A [TeV] 



Fig. 5 Expected and observed 95 % CL lower limits on the 
minimal GMSB model parameters A and tan /3. The dark 
grey area indicates the region which is theoretically excluded 
due to unphysical sparticle mass values. The different NLSP 
regions are indicated. In the CoNLSP region the ri and the 
£r are the NLSPs. Additional model parameters are Mmess 
= 250 TeV, Ns = 3, fi> and Cgrav = 1- The limits from the 
OPAL experiment [25] are shown for comparison. The recent 
ATLAS limit [22] obtained on a subset (2 fb-^) of the 2011 
data in the 2t final state is also shown. 



In the four final states studied, no significant excess 
is found above the expected SM backgrounds. The re- 
sults are used to set model-independent 95 % CL upper 
limits on the number of signal events from new phe- 
nomena and corresponding upper limits on the visible 
cross-section for the four different final states. Limits 
on the model parameters are set for a minimal GMSB 
model. A lower limit on the SUSY breaking scale A 
of 54 TeV is determined in the regions where the f i 
is the next-to-lightest SUSY particle (tan/3 > 20) by 
statistically combining the result of the four analyses 
described in this paper. The limit on A increases to 
58 TeV for tan (3 between 45 and 55 . These results pro- 
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vide the most stringent test to date of GMSB SUSY 
breaking models in a large part of the parameter space 
considered. 
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S.S. Bocchetta™, A. Bocci^^^ C.R. Boddy"^ M. Boehler^^, J. Boek^^^^ N. Boelaert3^ J.A. Bogaerts^o, 
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K. Cranmerio*, F. Cresciolii22a,i22b^ M Cristinziani^i, G. Crosetti^^^^'^^b^ s. Crepe-Renaudin^s, C.-M. Cuciuc^^'^, 

C. Cuenea Almenar^^^ T. Cuhadar Donszelmann"9, M. Curatolo^^, C.J. Curtis^^, C. Cuthbert^™, 
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R. Di Sipio20'*'20b^ m.A. Diaz32a, E.B. Diehl^^, J. Dietrich42, T.A. Dietzsch^^*^, S. Diglio*^ K. Dindar Yagci-^o, 
J. Dingfelder2i, F. Dinut26^, C. Dionisii32a,i32b^ p p)ita26a^ g Dita^S'^, F. Dittus^", F. Djama83, T. Djobava^i'', 
M.A.B. do Vale24<=, A. Do Valle Wemansi24a,n^ T.K.O. Doan^, M. DobbsS^ D. Dobos^o, E. Dobson^O-", 

J. Dodd3^ C. Doglioni^^, T. Doherty^^^ Y. Doi^^'*, J. Dolejsi^^e^ j Dolenc^^, Z. Dolezal^^e^ g a. Dolgoshein^^^-*, 
T. Dohmae^ss^ M. Donadelli^^d, j. Donini^^, J. Dopke^^, A. Doria^^^a^ a. Dos AIljosl^^ A. Dottii22a,i22b^ 



14 



M.T. Dova™, A.D. Doxiadisi"'\ A.T. Doyle^^ N. Drcssnandti^o^ M. Drisi", J. Diibbcrt^s, S. Dubci'\ 

E. Duchovni^'^^ G. Duckeck9^ D. Dude^'^^, A. Dudarev^^, F. Dudziak'^^^ jyj Dlihrssen^o, LP. Duerdoth*^^ 

L. Duflot"^, M-A. Dufour^s^ l Duguid^^ M. Dunford^^'^, H. Duran Yildiz^'^, R. Duxfield^^'', M. Dwuznik^^, 

F. Dydak^o, M. Diiren^^, W.L. Ebcnstcin^s, J. Ebkc'^^ S. Eckwcilcr^i, K. Edmonds^i, W. Edson^, 
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